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Abstract. Let A be a positive number, and let (xj : j £ Z) C R be a fixed 
Riesz-basis sequence, namely, (xj) is strictly increasing, and the set of functions 
{R 3in e lx i t : j £ Z} is a Riesz basis (i.e., unconditional basis) for La[— jr, 7r]. 
Given a function / 6 whose Fourier transform is zero almost everywhere 

outside the interval [— w, tt], there is a unique square-summable sequence (a,j : 
j £ Z), depending on A and /, such that the function 

i x (f)(x) := Y; a 3 e ~ Mx ~ Xj)2 > xeR, 

is continuous and square integrable on (—00,00), and satisfies the interpolatory 
conditions I x {f)(xj) = f(xj), j £ Z. It is shown that I x (f) converges to / 
in L2QR), and also uniformly on R, as A — > + . A multidimensional version 
of this result is also obtained. In addition, the fundamental functions for the 
univariate interpolation process are defined, and some of their basic properties, 
including their exponential decay for large argument, are established. It is 
further shown that the associated interpolation operators are bounded on £ P (Z) 
for every p £ [1, 00]. 

1. Introduction 

This paper, one in the long tradition of those involving the interpolatory theory 
of functions, is concerned with interpolation of data via the translates of a Gaussian 
kernel. The motivation for this work is twofold. The first is the theory of Cardinal 
Interpolation, which deals with the interpolation of data prescribed at the integer 
lattice, by means of the integer shifts of a single function. This subject has a rather 
long history, and it enjoys interesting connections with other branches of pure and 
applied mathematics, e.g. Toeplitz matrices, Function Theory Harmonic Analy- 
sis, Sampling Theory. When the underlying function (whose shifts form the basis 
for interpolation) is taken to be the so-called Cardinal B-Spline, one deals with 
Cardinal Spline Interpolation, a subject championed by Schoenberg, and taken up 
in earnest by a host of followers. More recently, it was discovered that there is a 
remarkable analogy between cardinal spline interpolation and cardinal interpola- 
tion by means of the (integer) shifts of a Gaussian, a survey of which topic may 
be found in |RS3j . The current article may also be viewed as a contribution in 
this vein; it too explores further connections between the interpolatory theory of 
splines and that of the Gauss kernel, but does so in the context of interpolation 
at point sets which are more general than the integer lattice. This brings us to 
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the second, and principal, motivating influence for our work, namely the researches 
of Lyubarskii and Madych [LMj . This duo have considered spline interpolation 
at certain (infinite) sets of points which are generalizations of the integer lattice, 
and we were prompted by their work to ask if the analogy between splines and 
Gaussians, very much in evidence in the context of cardinal interpolation, persists 
in this 'nonuniform' setting also. Our paper seeks to show that this is indeed the 
case. The influence of [LM] on our work goes further. Besides providing us with 
the motivating question for our studies, it also offered us an array of basic tools 
which we have modified and adapted. 

We shall supply more particulars - of a technical nature - concerning the present 
paper later in this introductory section, soon after we finish discussing some requi- 
site general material. 

A basic tool in our analysis is the Fourier Transform, so we assemble some basic 
and relevant facts about it here; our sources for this material are |Go| and |Ch| . If 
g € Li(W), then the Fourier transform of g, g, is defined as follows: 

/oo 
g{t)e~ ixt dt, iel. 
-oo 

It is known that g is uniformly continuous on R, and that lmi x —>± 00 g(x) = 0. In 
general g need not be integrable, but if it is, and if g G C(R) (the space of functions 
which are continuous throughout the real line), then one obtains the following 
inversion formula: 



1 f°° 

9® = — I 9{x)e lxt dx, te 
2tt J^oo 



Suppose now that g £ L 2 (R). The Fourier transform of g, denoted by T\g], is 
the function in -L 2 (R) for which limAr^oo \\1in — J~[g] ||i 2 (R) = 0, where 



JV 

ixt 



h N (x) := / g{t)e lxt dt, xGR. 

J-N 

The integral above is finite for every real number x, because g is square integrable 
on R, hence locally integrable on R. As ^[g] is obtained (a priori) only as an 
element in L 2 (R), it is determined only almost everywhere. It is known that T is 
a linear isomorphism on L 2 (R), and that the following hold: 

(2) irailLm - 27r||<7||! 2(R) , g € L 2 (R); T[g\ =g,g£ L 2 (R) n L ± (R). 

The inversion formula for the Fourier transform of square-integrable functions takes 
the following form: limN^a© ||<? — #jv||z, 3 (k) = 0, where 

1 f N 

(3) H N (t):=— T[g]{x)e lxt dt, t G R. 

2tt J_ n 

If, in addition to being square integrable, T[g\ is also integrable on R, then the 
Dominated Convergence Theorem implies that 

1 1"°° 

(4) lim Hjf(t) = — / T[g]{x)e lxt dx =: H(t), iel, 
n^oo 2tt J_ 00 
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so g must coincide with H almost everywhere. As H is continuous on R (in fact, 
2-KH{t) = F\g] (-t)), we find that, if g£L 2 (M.) n C(R) and T[g\ G£i(R), then 

(5) g{t) = 2^J ?\9K x ) e **> teM - 

The functions we seek to interpolate are the so-called bandlimited or Paley- 
Wiener functions. Specifically, we define 

PW T := {<? € £2(K0 : J~[g] = almost everywhere outside [— ir, 7r]}. 

The first equation in |J2) leads to the finding that PW V is a closed subspace of L2 (R) . 
Moreover, if g G PW^, then T[g] G .^(R) and .F[</] = almost everywhere outside 
[— 7r, 7r], so J-[g] G Li(R); hence the inversion formula discussed in the foregoing 
paragraph asserts that 



(6) g® = I 



T[g] {x)e lxt dx = —J T[g] (x)e lxt dx, 



for almost every real number t. Therefore, by altering the values of g, if need be, 
on a Lebesgue-null set, we may assume that ([6]) holds for every real number t, 
and this assumption will be in place throughout the article. In particular we shall 
assume that every function in PW n is continuous throughout R. Moreover, the 
Bunyakovskii-Cauchy-Schwarz Inequality and ^ combine to show that g is also a 
bounded function: 

(7) l5(*)l<-iiraiU a [-^] = NI^(R), teR. 



Even though it is not relevant here, we mention, at least by way of explaining our 
terminology, that the celebrated Paley- Wiener Theorem proclaims that a function 
belongs to PW^ if and only if it can be extended to the complex plane as an entire 
function of exponential type at most it. 

Having discussed the functions which are to be interpolated, we now describe the 
(canonical) point sets at which these functions will be interpolated; it is customary 
to refer to these points as interpolation points, or sampling points, or data sites. For 
the most part, though not always, we shall be concerned with data sites which give 
rise to Riesz-basis sequences. Precisely, following [LMJ, we say that a real sequence 
(xj : j G Z) is an Riesz-basis sequence if it satisfies the following conditions: Xj < 
xj + i for every integer j, and the sequence of functions (ej(t) := e~ lXjt : j G Z, t G 
R) is a Riesz basis for i 2 [~ 7r : ""]• We recall that saying that a sequence (tpj : j G Z) 
in a Hilbert space 7i is a Riesz basis for Ti. means that every element h G TL admits 
a unique representation of the form 

(8) h = J2 a 3<fji EM 2 <oo, 

and that there exists a universal constant B such that 

1/2 



o) b- 1 hr 



< 



j'ez 



< £ 



H 




for every square-summable sequence (cj : j G Z). Classical examples of Riesz-basis 
sequences are given in |LM| . where it is also pointed out that if (xj : j G Z) is an 



4 



TH. SCHLUMPRECHT AND N. SIVAKUMAR 



Riesz-basis sequence, then there exist positive numbers q and Q such that 

(10) q < Xj + i — Xj < Q, j £ Z. 

The interpolation process we study here is one that arises from translating a 
fixed Gaussian. Specifically, let A > be fixed, and let (xj : j £ Z) be an Riesz- 
basis sequence (a relaxation of this condition will be considered in the final section 
of the paper). We show in the next section that given a function / £ PW^, there 
exists a unique square-summable sequence (<ij : j £ Z) - depending on A, /, and 
the sampling points (xj) - such that the function 

(11) I x (f)(x):=Y,a 3 e-^-^\ x £ R, 

jez 

is continuous and square-integrable on R, and satisfies the interpolatory conditions 

(12) I x {f)(x k ) = f{x k ), fceZ. 

The function I x (/) is called the Gaussian Interpolant to / at the data sites (x^ : 
k £ Z). We also prove, again in the upcoming section, that the map / h- > I x (f) 
is a bounded linear operator from PW V to £2 (R)- As expected, the norm of this 
operator I x - which we refer to as the Gaussian Interpolation Operator - is shown 
to be bounded by a constant depending on A and the choice of the Riesz-basis 
sequence. However, this is not sufficient for our subsequent analysis, in which we 
intend to vary the scaling parameter A. So in Section 3 we demonstrate that, if the 
underlying RRB sequence is fixed, and if A < 1 (the upper bound 1 being purely 
a matter of convenience), then the operator norm of I x can be majorized by a 
number which is independent of A. Armed with this finding, we proceed to Section 
4, wherein we establish the following focal convergence result: 

Theorem 1.1. Suppose that (xj : j £ Z) is a (fixed) Riesz-basis sequence, and let 
I x be the associated Gaussian Interpolation Operator. Then for any f £ PW^, we 
have f = lim^^g+ I x {f) in 2v2(R) and uniformly on R. 

We note that, in the case when Xj — j, this theorem was proved in |BSj . 

Our proofs in Sections 2-5 rely heavily on the machinery and methods developed 
in |LMj for cardinal splines; indeed, as mentioned earlier, our primary task in this 
paper has been to adapt these to the study of the Gaussian. However, most of 
these arguments do not extend per se to the multidimensional situation, which 
occupies our attention in Section 5. The results presented in this section are far from 
complete, and should be viewed only as partial generalizations of their univariate 
counterparts. Nonetheless, it is not without interest to note that tackling even 
this simplified situation requires a combination of the results established in one 
dimension and some abstract functional-analytic techniques. The paper concludes 
with Section 6, in which we revisit univariate interpolation, but consider sampling 
points which satisfy a less restrictive condition than that of giving rise to a Riesz- 
basis sequence. We introduce here the fundamental functions for interpolation at 
such data sites, and prove that they decay exponentially for large argument. In 
addition to being of independent interest, as readers familiar with spline theory 
will readily attest, this result also paves the way towards a generalization of some 
of the main results of Section 2. 

We have attempted to make this article as self contained as possible, and we re- 
quest the indulgence of those readers who may find an abundance of detail between 
these pages. 
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2. NOTATIONS AND BASIC FACTS 

In this section we shall reintroduce the interpolation problem which concerns 
us, define the corresponding interpolant and interpolation operator, and establish 
some of their basic properties. We shall uncover these in a series of propositions, 
which begins with this simple observation. 

Proposition 2.1. If a is a positive number, then 

E e- Q(2 " hl)2 < ^ =: «(«). 
zez\{o} 

Proof. 

^2 e -«(2i/i-D 2 < 2 jr e- al = 

iei,\{o} i=i 6 

□ 

In what follows we shall use the following notation: given a positive number A, 
the Gaussian function with scaling parameter A is defined by 

g x (x) := e- Xx \ i£l 
We recall the well-known fact (see, for example, |Go[ p. 43]) that 

(13) ^ A ](«) = Tx{u) = y| e -« 2 /( 4A ), «ei. 

We now record two results from the literature; both will be of use in this section 
and also in Section 5. 

Proposition 2.2. cf. [NSW1 Lemma 2.1] 

Let A and q be fixed positive numbers, and let \\ ■ \\ 2 denote the Euclidean norm in 
Mr. There exists a number v, depending only on d, X, and q, such that the following 
holds: if (xj) is any sequence in M. d with \\xj — Xk\\ 2 > q for j ^ k, and x is any 
point in M. d , then g\ ( 1 1 x — Xj\\ 2 ) < v. 

This next result is an important finding in the theory of radial-basis functions. 

Theorem 2.3. cf. [NWl Theorem 2.3] 

Let A and q be fixed positive numbers, and let || • ||a denote the Euclidean norm 
in M. d . There exists a number 9, depending only on d, X, and q, such that the 
following holds: if (xj) is any sequence in M. d with \\xj — Xk\\2 > q for j 7^ k, then 
J2j,k^k9\(\\xj - x k \\ 2 ) > 9J2j ICj| 2 ; for every sequence of complex numbers 

Proposition 2.4. Suppose that (xj : j G Z) is a sequence of real numbers satisfying 
the following condition: there exists a positive number q such that > a 

for every integer j . Let A > be fixed, and let (a,j : j € Z) be a bounded sequence 
of complex numbers. Then the function R 3 x 1— * X^/ez a j9^( x ~ x j) * s continuous 
and bounded throughout the real line. 



Proof. Proposition 12.21 demonstrates that the series in question is uniformly con- 
vergent throughout R, and that the sum is a bounded function of x. The apparent 
continuity of each summand and uniform convergence imply the continuity of the 
limit function. □ 
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Remark 2.5. Suppose that A is a fixed positive number, and let (xj : j G Z) 
be a sequence satisfying the conditions of Proposition 12.41 The latter conditions 
imply that \xj — x^\ > \j — k\q for every pair of integers j and k, so the entries 
of the bi-infinite matrix (g\(xk — Xj))k.jez decay exponentially away from its main 
diagonal. So the matrix is realizable as the sum of a uniformly convergent series of 
diagonal matrices. Hence it acts as a bounded operator on every ^ P (Z), 1 < p < oo. 
Moreover, as the matrix is also symmetric, Theorem l2.3l ensures that it is boundedly 
invertible on £2 (20- In particular, given a square-summablc sequence (dk ■ k G Z), 
there exists a unique square-summable sequence (a(J, A) : j G Z) such that 



X)g\(x k - Xj) = d k , 



fc g z. 



Suppose now that (xj : j £ Z) is a Riesz-basis sequence. Thus, given h G 
L2[— 7r, 7r], there exists a square-summable sequence (aj : j g Z) such that /i(t) = 



cue 



as follows: 



for almost every i G [— 7T, 7r]. We wish to extend this function to 



Proposition 2.6. Let (xj : j G Z), /i, and (oj : j £ Z) 6e as above. The function 
H(u) := J2j£Z a j e ~ iXjU * s l oca Hy square integrable on M; m particular it is well 
defined for almost every real number u. 



Proof. Recall from the introductory section that there exists a constant B such 
that 



(14) 



jez 



jez 



dt<B 2 J2\ 



for every square-summable sequence (cj : j 6 Z). Let u G [(2/ — l)7r, (2Z + l)7r] for 
some integer I, and define, for every positive integer N, the (continuous) function 



H N (u) := E^-jv a i< 

f (2i+l)7T 



(15) 



If N > M are positive integers, then 

2 



(2Z-1)tt 



(2Z+1)tt 
(2;-1)tt 



AT 

E °j' 

|j|=Af+l 



iV 



< £ 



b'|=M+i 

JV 

E l< 

b'|=M+i 



a 7 e J e J 



N 

E 

|j|=M+l 



the final inequality coming from f| 14[) . Ergo, the square summability of the sequence 
(aj : j G Z) shows that (iijv : G N) is a Cauchy sequence, hence convergent, in 
L2[(2l — l)-7r, (2Z + l)7r]. This proves the desired result. □ 

For future reference, we note that the argument leading up to (|15p also provides 
the estimate 



(16) 



n (2; + l)ir 



\H\ 



'(2Z-1)7T 



/•(2/+l)jr 

lim / 

JV - >00 7(2/-l)7T 



Iff 



< 



z G z, 
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whence the Bunyakovskii-Cauchy-Schwarz Inequality implies that 

1/2 



(17) 



,(2Z + 1)tT , 

/ \h\<V2^b\J2\ 

J(2l-l)ir \ jez 



I G Z. 



The next result is the first of the two main offerings of the current section. 

Theorem 2.7. Suppose that X is a fixed positive number, and let (xj : j 6 Z) be a 
Riesz-basis sequence. Assume that a :— (aj : j G Z) is a square- summable sequence. 
The following hold: 

(i) The function 

s(a,x) = s(x) :='^^ajg\(x — Xj), x G M, 
j&L 

belongs to C(E) H L 2 (M). 

(ii) The function 

S(a,u) = §{u) := e- u2 /^J2 a J e ~ ZX3U 
is well defined for almost every real number u, and s G L2(R) H Li(R). 



is a bounded 



(iv) The map a :— (aj : j £ Z) h s(a,ir) := J2jez a j9><( x ~ x j)> x ^ 
linear transformation from £ 2 (Z) into L2(R). 



Proof, (i) The continuity of s on K follows at once from the first inequality in |T 



and Proposition 12.41 Define sn(x) :— J2j=-N a j9^( x ~ x j)> i € R, JV £ N. As 
sn G Li(M)ni2(K) for every positive integer N, equation (flU)) . and a standard 
Fourier-transform calculation provide the following relations: if N > M are positive 
integers, then 



(18) 2ir\\s N - s M ||L 



II^V ~ •H/||| 2 (R) 

p -« 2 /(2A) 



E • 

b'l=M+i 



(2i+l)7T 



-« 2 /(2A) 



; £Z J(22-1),T 



AT 

E 

\j\=M+l 



a^e 



du. 



Let H m ,n{u) := J2 



N 

j\=M+l 



G R. Using the estimates e"" 2/(2A) < 1 



for u G [-7T,7r], and e~ u < e — i| — i) -k /(2A) for ^ _ 1)^ < u < ( 2 / + 1)tt, 
Z G Z \ {0}, we find from OH) that 



(19) 2ir\\s N -s M 



< 



77 r ,, 

- X [\\H M 



N 



L 2 [— 7T,7r] 



4- V e -(2|'|-l) 2 ^/(2A) ||rr ||2 

;ez\{o} 
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Now let B be the constant satisfying (|14|) and use (fT5|) to estimate each of the 
integrals on the right side of (|19|) . This yields the relation 



27r||sAr - s M \\ 



< 



B 2 



N 

E 



2 



Jj|=M+l 

whence Proposition 1 2 . 1 1 leads to the estimate 

rl- ^ / N 



E < 

zez\{o} 



-(2|;|-1)V 2 /(2A) 



27r||s7v - s M \\ L2 ( R) < 



B IT 



E Kl 2 ] [1 + ^V(2A))] 

|j|=Af+l 



As (oj : j £ Z) is square summable, we find that (sn ■ N £ N) is a Cauchy sequence 
in La(K), and hence that s € as promised. 

(ii) Let H(u) := /Cjez a,je~ lXjU . As observed in Proposition 12.61 H is defined 
almost everywhere on K, so the same is true of s as well. Now the argument in (i), 
combined with (TIB)) , shows that 



(20) 



~|| 2 



< B2 EN' [! + ^7(2A))] 5 



whilst a slight, but obvious, variation on the theme, coupled with (|17|) . demonstrates 
that 



1/2 



(21) 



|5|U lW < V^B Yl Kf [1 + ^V(2A))] , 



and this completes the proof. 

(in) Let (sjy : AT e N) be the sequence defined in the proof of (i). As each 
sn € n an d limjv-*cw H s iV — s IIl 2 (R) = 0, ^ suffices to show, thanks 

to @, that limjv^oo || Sjv — s/ (tt/A) s||l 2 (r) = 0. Calculations similar to the one 
carried out in (i) show that 



SAT 



-ii 2 



-« 2 /(2A) 



E 



b1=iv+i 

2 



du 



< 



B 2 TT 



A 



E MM [l + «(^/(2A))] 



and the last term approaches zero as N tends to infinity, because the sequence 
(aj : j £ Z) is square summable. 

(iv) The linearity of the map is evident, and that it takes £2(2) into is 
the content of part (i). Now ((2]), part (iii) above, and (f20|) combine to yield the 
relations 



27r|| S (a,.)|| 



J s(a,-) 



< 



B 2 TT 



1 + K (^ 2 /(2A))] ||ft 



2 

1*2(2) 



□ 
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This next result points to a useful interplay between Riesz-basis sequences and 
bandlimited functions (see, for example, |Yo[ pp. 29-32]). It serves as a prelude to 
the second main theorem of this section. 

Proposition 2.8. Suppose that (xj : j G Z) is a Riesz-basis sequence, and that 
f G PWtt. Then the (sampled) sequence (f(xj) : j G Z) is square summable. 
Moreover, there is an absolute constant C - depending only on (xj), but not on f 
-such that E 7 «I/(^)| 2 <C 2 ||/||2 . 



Proof. Let 
(22) 



(hi, hi 



hih 2 



hi, h 2 e L 2 [-7r,7r] 



denote the standard inner product in L 2 \— tt, tt]. Let ej(t) 



f , j ez,te 



[— 7r,7r], so that ([6]) implies the identities 

(23) 2TTf(x j ) = (F[f],e j ), jeZ. 

Letting (ej : j G Z) be the co-ordinate functionals of (e^ : j £ Z) (which means 
that h — Y]j(h, e,-)e^ for any h G L 2 [7r,7r]), it follows that (e.j : j 6 Z) is also a 
Riesz basis whose co-ordinate functionals are (e^ : j 6 Z). Thus, 

(24) g = ^2(g,ej)ej, g G i 2 [-tt, tt] . 



So ([9]) provides a universal constant -B such that 



Ei 



< £r 



E 

j6Z 



i2[ — 7T,7r] 



for every square summable sequence (cj : j G Z) . Hence (|2"3"|) and (f2~4"|) imply that 

4^ 2 ^ i/(^-)i 2 < wmwU-*,*] = 2^ 2 n/ni 2(R) , 

the final equation stemming from ([2]), and the fact that T[f] — almost everywhere 
in R\ [-7T,7r]. □ 

We now state the second of the two main results in this section. Most of our 
work has already been accomplished; what remains is to recast the findings in the 
context of our interpolation problem. 

Theorem 2.9. Let A be a fixed positive number, and let (xj : j G Z) be a Riesz-basis 
sequence. The following hold: 

(i) Given f G PW^, there exists a unique square- summable sequence (a(j,\) : j G 
Z) such that 

E a ^' X )9*( x k ~ Xj) = f(x k ), keZ. 
jez 

(ii) Let f and (a(j, A) : j G Z) &e as m (%). T/ie Gaussian interpolant to f at the 
points (xj : j G Z), to wif, 

h(f)( x ) =z2a(j,><)9x(x - xj), x G M, 
jez 

fceZor^s to C(K)nL 2 (l). 
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(Hi) Let f and I x {f) be as above. The Fourier transform of I x {f) is given by 



F\h (/)](«) = 




for almost every real number u. Moreover, ZF[I x (f)] G H Li(R). 

(zt;,) // / and I x {f) are as above, then 

hif)(x) = f°° F[i x (f)](uy™du, xeR. 
*K J -00 

In particular, I x (f) £ C (R) := {g G C(M) : lim^i^oo = 0}. 
(v) The Gaussian interpolation operator I A is a bounded linear operator from PW^ 
to L/2(M.). That is, the map PW^ 9 / 1— > I x (f) is linear, and there exists a positive 
constant D, depending only on A and (xj : j G Z), such that 

II^a(/)IIl 2 (e) < -D II/IIl 2 (r) 

for every f £ PW^ . 

Proof. Assertion (i) follows from the first inequality in (fTO]) . Remark 12.51 and 
Proposition 12.81 Assertion (ii) obtains from part (i) of Theorem 12. 7\ whilst as- 
sertion (iii) is a consequence of parts (ii) and (hi) of Theorem 12.71 Assertion 
(iv) follows from the fact that I x (/) is continuous throughout R, that J-[I X (/)] € 
Z/2(R) n .Li(R), and equation ([5]). Moreover, this representation and equation ((T|) 

show that 2irl x (f)(x) = J-[I x (f)](—x) for every real number x. So the Riemann- 
Lebesgue Lemma |Go[ Theorem 4A] ensures that I x € Co(R). As to (v), let T be 
the matrix (g\(xk — Xj))j t kez, let / £ PW^, and let d :— (f{xj) : j £ Z). Then 
a(j, A) is the j-th component of the vector T~ 1 d, and this demonstrates that I x is 
linear. Now part (iv) of Theorem 12 . 71 asserts that 

(25) IK(/)II M m) = O (IIT-^I^^) = O (|R||, 2(Z) ) , 

where the Big-0 constant depends only on A and the Riesz-basis sequence (Xj : j £ 
%). Furthermore, Proposition 12 . 81 reveals that 

( 26 ) IRL(z) = (ii/IIl 2 ( R) ), 

with the Big-0 constant here depending only on (xj : j G Z). Combining (|25p with 
(|26| finishes the proof. □ 

3. Uniform boundedness of the interpolation operators 

In the final result of the previous section, it was shown that, for a fixed scaling 
parameter A, and a fixed Riesz-basis sequence (xj : j £ Z), the associated interpo- 
lation operator I x is a continuous linear map from PW„ into La(R). As expected, 
the norm of this operator was shown to be bounded by a number which depends 
on both the scaling parameter and the choice of the Riesz-basis sequence. The goal 
in the current section is to demonstrate that, if the scaling parameter is bounded 
above by a fixed number (taken here to be 1 for convenience) , then the norm of I x 
can be bounded by a number which depends only on {xj : j G Z). The proofs in 
this section (as well as in the next) are patterned after |LM) . 

Let (xj : j £ Z) be a Riesz-basis sequence, and let B be the associated constant 
satisfying the inequalities in (fl4|) . Given h £ L^y—it, tt], there is a square-summable 
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sequence (<Zj : j G Z) such that h(t) = J2jez a j e ~ lXjt f° r almost every t G [— 7r, 7r]. 
Let H denote the extension of h to almost all of K, as considered in Proposition 
Given an integer I, we define the following linear map Ai on L 2 [— tt,tt]: 



(27) Ai (h) (t) := H(t + 2nl) = ^ aj e^ (t+27ri) 

for almost every t G [— tt,tt]. We see from p^|) and (jT^J) that 

(28) P/^llL^j = l|ff|ll [( 2/-iK(2 i+ i).] < B2 E M' ^ 54 ll^llLh^] ■ 

Thus every Ai is a bounded operator from Ln\— 7r,7r] into itself; moreover, the 
associated operator norms of these operators are uniformly bounded: 

(29) \\M\\<B 2 . 

In what follows, we shall assume that (xj : j G Z) is a (fixed) Riesz-basis 
sequence, and let Cj(t) := e _M:jt , { £ 1, j 6 Z. We also denote by (•,•) the 
standard inner product in L^y—i^ as defined via (|22|) . Our first main task now 
is to exploit the presence of the Riesz basis (ej : j G Z) in L2[— 7r,7r] to find an 
effective representation for the Fourier transform of the Gaussian interpolant to 
a given bandlimited function, on the interval [— 7r, tt]. We begin with a pair of 
preliminary observations: 

Lemma 3.1. Let (xj) and (ej) be as above, and let f be a given function in PW n . 
If <p G Li2[— tt, tt] satisfies the conditions 

/IT 
c/)(t)e~ iXkt dt, fcGZ, 
-7T 

i/ien ■?""[/] agrees with <f> in L%\— 7T, tt]. 

Proof. Equations (|23p and ([30|) reveal that (^F[/],efc) = ((f), ek) for every integer fc, 
and the required result follows from □ 

Lemma 3.2. Lei (xj) and (ej) be as above, and let B be the constant satisfying 
(fH| . Let h G L-2 [— 7T, tt] , and Zet a > 0. Define 



<j>l=A]{e-^ + ^A l (h)), 

Then 



fcGZ. 



(31) ||0o|U 2 Hr,7r] < IHU 2 [-7r,7r] «^ E H ^ H £2 ^ II fe H L 2 [-tt.tt] B i n(lT 2 a) , 

;ez\{o} 



where k is the familiar function from Proposition \2.1\ In particular, the series 
^2i e %4>i converges in Li\— 7T, tt]. 

Proof. We note that 

(32) lltoll^l-*,,] = < II^IU,^,,], 
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whereas (|29| gives rise to the following estimates for every i 6 Z \ {0}: 



(33) 



''HW-TT.TI-l 



<B 1 



Lo\-T. 



<B*e-^ l \-V 2 \\h\\ L ^ M . 

This completes the proof. 

We are now ready for the first main result of this section. 



□ 



Theorem 3.3. Let A > be fixed, and let f S PW^. Let ip\ denote the restriction, 
to the interval [— it, ir], of the function "J a given in part (Hi) of Theorera \2.9\ Then 

(34) m=FMf)] + M E Al^+W^Aityxj) on [-*,*]. 

" iez\{o} 

Proof. Let <j> denote the function on the right side of (f3"4"|) . In view of Lemma 13.11 
it suffices to show that 

(35) 2nf(x k ) = {<j>, e k ), k e Z. 
Now Theorem [279] implies the relations 

(36) 2nf{x k ) = 2ttI x (f)(x k ) = T e~ u * ' ^ x {u)e™^ du, k e Z, 



whilst 



(37) 



,-« 2 /(4A) 



/•(2i+l)ir 

^>x{u)e lXkU du = J2 e- u ^ 4X ^ x (u)e lXkU du 

= J2 r e -(™) 2 /( 4A )* A (t + 2wl)e ix ^ t+M Ut 
lei , n 

= T,r ^ it+M)2/(4X) MM(t)A^W)dt 

E(A*(e- ( - +M)2/(4A) ^(^)),e fe ; 
Y^A] {e-^^I^M^)) ,e k 



the final step being justified by Lemma [3. 21 Noting that 

J2 A "i (e-^ 2 /^AM)=e-^ 2 ^ x + £ A* (e-(-+^) 2 /( 4A )A ; (^)) 



/e 



iez\{o} 

we find that [[3T]). ([315]). and part (iii) of Theorem [231 yield (J35J), and with it the 
proof. □ 

Combining ([3"4"|) and ([3"Tj) leads directly to the following: 
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Corollary 3.4. Suppose that n, A, /, ip\, and B are as before. Then 

II^(/)]Hl 2[ -^] < WmWL,^ + \[\b^{tt 2 /{±\)) UxW^-^ ■ 

The preceding corollary shows that, if / is bandlimited, then the energy of the 
Fourier transform of (its Gaussian interpolant) I x (/) - on the interval [— 7r, tt] - is 
controlled by that of the Fourier transform of / on that interval, plus another term 
which involves the energy of ip\ on the interval. Our next task is to show that 
this second term can also be bounded effectively via the energy of T[f] on [— tt, tt]. 
Before proceeding with this, however, we pause to consider formally, an operator 
which has already made its debut, albeit indirectly, in Theorem 13.31 This operator 
will also play a role later in this section, and a larger one in the next. 

Given a positive number a, we define the operator T., on L2[— tt, tt] as follows: 

(38) T [a] (h):=e^ a ]T A\ (e^^'J^W) , h G L 2 [-tt, tt] . 
iez\{o} 

That this operator is well defined is guaranteed by Lemma T3.21 whilst its linearity 
is plain. The following properties of T,, are easy to verify. 

Proposition 3.5. The operator T.. is self adjoint, positive, and its norm is no 

larger than e v a B 4 K(TT 2 a), where k is the function defined through Proposition ^. 1\ 
and B is the familiar constant associated to the given Riesz-basis sequence (xj). 

We now return to the task of carrying forward the estimate in Corollary 13.41 
The first order of business is to attend to ||^A|li 2 [_ 7r 

Proposition 3.6. The following holds: 

\m^ { -^<yll^ m) \\m\\L^M- 



Proof. Equation (|34|) asserts that 

(39) m =F[lAf)\ +^ 7r2/(4A) /f ^/(^(V'a). 
As 

(40) OT*(/)Ua) = {{ff u2/i4X) Mu) ? du ^ °< 

and T [1/{ix)] is a positive operator, we find from (|39jl that {J-[f],%p\) is nonnega- 
tive, and also that (3~\f\,ip\) > (J~\I X (/)], "4>x)- Hence the Bunyakovskii-Cauchy- 
Schwarz inequality and (|40|) lead to the relations 

\\m\\M-*M \Mi 2 \-^ > \/f f e _t,V(4A) \Mu)\ 2 du 



>\^ 2/m \\Ml 2l ^ v 



A 

and the required result follows directly. □ 

The upcoming corollary is obtained via a combination of Corollary 13. 4[ Propo- 
sition[321 Equation ^ and the fact that e^ ' ' ' 4A k(tt 2 / (4A)) = (1 - ! /4A)-i < 2) 
whenever A < 1. 
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Corollary 3.7. Assume that < A < 1, and let (xj) and f be as above. The 
following holds: 



mh(f)}\\ L2l -^ ] <v^[i + 2B*]\\f\\ 



L-2 



The preceding result accomplishes the first half of what we set about to do in 
this section. The second part will be dealt with next; our deliberations will be quite 
brief, for the proof is now familiar terrain. 



Proposition 3.8. Let < A < 1. The following holds 



nhU)\\\L^\[-,M) ^ B\r-U x \\ L2[ _^n{^/{2\)) < V^B 2 \\f\\ 



Proof. We begin by noting that 

(41) WFMMl^-^) = \l e-» 2 /^\V x (u)\ 2 du 

/ .(2/+1)tt 

= ? £ / e- 2 /^ x (ufdu. 
The last term in (|4"Tj) may be bounded as follows: 



( 42 ) I £ ff^ 2 /^\(AM,))(t)\ 2 dt 

<- x E e - (2|;hl)V2/(2A) p i (^)ii 2 L2[ -^ ] 



< 



zez\{o} 
B 4 n 



MxWU-^] «(^/(2A)), 



A 

the last inequality being consequent upon (I29j) . This proves the first of the two 
stated inequalities. 

The second inequality follows from the first, by way of Proposition 13. 6[ (J2J) and 
the fact that (1 - e -^ 2 /4A-)-i < 2j whenever A < 1. □ 



We close this section by summarizing the findings of Corollary 13. 71 and Proposi- 
tion [37 



Theorem 3.9. Suppose that (xj : j S Z) is a fixed Riesz-basis sequence. Then 
{I x : < A < 1} is a uniformly-bounded family of linear operators from PW^ to 
L 2 (R). 

4. Convergence of I x 

This section is devoted to the proof of the convergence result stated in the 
introduction (Theorem 1 We begin by laying some requisite groundwork. Let 
a be a fixed positive number. Recall the linear operator T [a] from (|38|) : 

T [a] (h):=e* 2a E A 'l (e- a{ - +27rl)2 Mh)) , h € L 2 [-tt, tt]. 
iez\{o} 

We now define the following (multiplier) operator on L^y— 7r, tt]: 
(43) M [a] {h) := e'^'^h, h G L 2 [-tt,it]. 

The following properties of M„ are easy to verify. 
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Proposition 4.1. The operator -M,, is a bounded linear operator on L2[—tt,tt], 
whose norm does not exceed 1. Moreover, it is self adjoint, strictly positive, and 
invertible. 

In what follows, we let (xj : j S Z) be a fixed Riesz-basis sequence, and let 
r : L2(M) — > Li\— ir,ir] denote the map which sends a function in L2(M) to its 
restriction to the interval [— 7r,7r]; note that r is a bounded linear map with unit 
norm. 

Let / £ PWn- Recall (from part (iii) of Theorem l2.9[) that, if J A (/) denotes the 
Gaussian interpolant to / at the points (xj), then 

(44) Mt) = VNv)e t2/{AX) niAf)]{t) 

for almost every t in [— tt, tt] (remembering that ip\ is the restriction of ^\ to the 
interval [— tt, it]). With all this in mind, we find from the definitions of T. , and 
■M [a] that equation may be cast in the following form: 

(45) (l + T [1/(ix)] M [1/(4X)] )r(T[I x (f)])=r(T[f]), 

where X denotes the identity on L2[~tt, tt]. 
Suppose that g 6 L2[—tt,tt], and let 

\g(t) if t e [— 7r,7r]; 
if t 6l\ [-7r,7r]. 



Then there is an / S PW^ such that T[f] = g; in fact, we may take / to be the 
following: 



g{t)e ixt dt, x e 



Let I x (f) denote the interpolant to / at the points (xj), and define L\(g) := 
r(T[I x (/)]); in other words, L\ = r o To I x o T~ x . As the Fourier transform T is a 
linear isomorphism on L2(R), and the maps I x and r are linear and continuous, the 
map g i— ► L\(g) is a continuous linear operator on L2[—tt,tt]; moreover, equation 
(l45l) affirms that 



(46) (l + T ll/(4X)] M ll/{4X)] )L x (g) = g. 

This being true for every g € L2[— 7r, 7r], we deduce that the map I+T, x ,, lx ,, -M [1/(4X)] 
is surjective on L2[— 7r, 7r], and that L\ is a right inverse of . (4X) . A4 (1 ,, 4A)] . Wc 

now show that X + M.,, is, in fact, invertible on L^f— 7r, 7rl. 

Proposition 4.2. TTie map X + T^.^. .M (1/(4A)] is injective, hence invertible, on 
L^[— 7T,7r]. Moreover, there is a constant A, depending only on the sequence (xj), 
such that 

p+t^m^)- 1 < a > o<a<i. 

Proof. Suppose that (l + 1~ [1/(ix)] A4 (1/(4X)] ) g = for some g G L2[— tt, tt]. Then 

= (( I + 7 i 1 /(4,) 1 >' [1/( 4„ ] ) (ff).M W(i X),(ff)) 

= fa. • M [1/ (4A) 1 (.9)) + (^ l/( 4A)]^[l/(*X)] M [1/(4X,] (9)) 

where we have used the positivity of the operators 7j 1/(4A)] and M [1/{4X)] to obtain the 
two inequalities above. It follows that M. [1/{ix)] {g) = 0, and, as VW [1/(4A)] is strictly 
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positive, g must be zero. Hence X + 7j 1/(4A)] VW [1/(4A)] is injective, and therefore 
invertible. So (|46|) may now be stated as follows: 



1 + T M 

V ^ jI [l/(4A)]-' l [l/(4A)]/ 



1 + T M 



for A G (0, 1] 



Consequently, the uniform boundedness of 

obtains from recalling the equation L\ = r o T o J o JF _1 , along with Theorem 
EH □ 



We are now ready for the first of the two focal results of this section. 
Theorem 4.3. If f e PW^, then 

fon ||/-/ A (/)|| i2(R) =0. 

Proof. In view of the first identity in ([2]), it is sufficient to show that 



(47) 



lim 



L 2 (R) 



= 0. 



Assume that < A < 1, and let A' := 1/(4A). As T[f] is zero almost everywhere 
outside [— 7r, 7r], we see that 

(48) \\T[f] -T[I x (f)]\\l 2m = \\m -^^(/)]l| 2 L 2[ - 7r ,. ] +ll^[^(/)]l| 2 L 2 ( R \ [ -., 7r] ) ■ 
On the interval [— tt,tt], we have, via (|43|) . that 



^[/]-^x(/)] = 



W] 

^^^^i^X^M WD, 

where the second step is a matter of direct verification. Consequently 
(49) ||^[/]-^[/ A (/)]|U 2 [-^] 



< 



WD 



L 2 [— 7r,7r] 



Now Proposition [472] provides a positive constant A which bounds the first term on 
the right of (ggj for every A e (0, 1]. As k(tt 2 /(4A)) = O^ 7 ^ 4 ^) for < A < 1, 
Proposition 13. 51 implies that 



O 



(e' r2 /( 4A) K (^ 2 /(4A))) = 0(1), < A < 1, 



for some Big-0 constant which is independent of A. Using this pair of estimates in 
(|49|) provides 



(50) \\m-F[hU)\\\ 



L 2 [— 7r,ir] 



= 



^,v, WD 



_L 2 [ — 7T ,7T 



< A < 1. 
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Turning to the second term on the right of (|48|) . we see from Proposition ^. 8( (|44|) . 
and (43|) that 



(51) ll^(/)]|lL ( E\[-^])=0 

= O 



lV'A|lL h ^ 1 ^( 7 r 2 /(2A)) \ 



L 2 [— tt.tt] 



K (rr 2 /(2A)) 



= e ^/(2X) 



= o 



M [XI] {F[I x (f)]) 



k (tt 2 /(2A)) 
< A < 1, 

2 - n^-^ 2 /(2A)\ 



the final step resulting from the (oft cited) estimate k(tt /(2A)) = 0(e 
< A < 1. Now 



(52) 



L 2 [— 7T,7r] 



M rvl (^[!x(/)]-nf]) 



L 2 f— 7r,7rl 



i 2 [— 7T,7r] 



L 2 [— 7r,7r] 



because < 1 (Proposition OJ, and ([50]) holds. Combining ([52]) and ijgQl 

with (JUJ, we find that 



(53) ||^[/]-^(/)]||L 



(n^ [1 /(.x )I w])ir i2[ -^])=°( 1 )' a 



0^ 



the last assertion being a consequence of the Dominated Convergence Theorem. 
This establishes (j47|) . and the proof is complete. □ 



The final theorem of the section deals with uniform convergence. 

Theorem 4.4. If f e PW^, then lim A ^ + h (f)( x ) = f( x ), x e R, and the 
convergence is uniform on R. In particular, the operators I\, < A < 1, are 
uniformly bounded as operators from PW n to Co(IR), via the Uniform Boundedness 
Principle. 



Proof. Assume that < A < 1, and let x € K. From ([6]), part (iv) of Theorem 1 2. 9 
and the fact that T[f] = almost everywhere on E \ [— n, it], we see that 

(54) f(x)-I x {f){x) = i- f -^[J A (/)]( u ))e» u d U 

27T./_^ 



1 

2^ 



^[J A (/)]( U )e iCT d W . 



The modulus of the first term on the right side of (|54p is no larger than 

(55) \\m-F[i,{f)\\\ Ll[ -. M =o( \\m-niAf)]h 2[ -,,,]) =*(i)> a-o+, 
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via (|2|) and Theorem 14.31 The second term on the right side of (|54|) is estimated in 
a familiar way: 



(56) 



T[I x {f)]{u)e lxu du 



< 



\T[I X (/)](«) | du 



T E 



(2J+1)tt 



- u2 /^\m x ( u )\du 



< 



iez\{o} J ( 2l - 1 ) n 
sj\ E e- (a|,| - 1)V/(4X) ll^(^)IL l[ _ iriir] 



zez\{o} 



0(7? E ^ (2 " hl)27r2/(4A) ii^(VA)ii L2[ _^ 

^A|| £2[ _^^(vr 2 /(4A)) ^ 



Borrowing the argument which led up to ([ST]) , (f5"2")l . and the final conclusion of (f5"3")l . 
we deduce that I\(f)(x) converges to f(x) uniformly in R. 

□ 



5. A Multidimensional Extension 

We consider now the multidimensional Gaussian interpolation operator. Let d £ 
N, and let (xj : j £ N) C R d . We say that (x 3 : j £ N) = £(j,2)> ■ ■ ■ , £(j,d) : 

j € N) C R d is a d- dimensional Riesz-basis sequence if the sequence (e") : j e N), 
with 

e : : [-7r,7r] d ^C, ejfa.ta, . . . ,t d ) ■= e~*<'^ = e"*^ »0.0*«, 

is a Riesz basis of L2[— Tj 7r] d . 

In general there is no natural indexing of the elements of {xj} by Z or Z d if 
d > 2; so we index generic d-dimensional Riesz-basis sequences by N. Later we 
shall concentrate on grids in R d , i.e., on Riesz bases in R d which are products of 
one-dimensional Riesz-basis sequences. In that case the natural index set is Z d . We 
note that, as in the 1-dimensional case, a Riesz basis-sequence in R d also has to be 
separated [LMj . 

The d- dimensional Gaussian function with scaling parameter A > is defined by 

g{ d) (x u x 2 , ...,x d ) = e- A M 2 = e ~ x ^U*i, x = ( Xl , X2 , . . . , x d ) £ R d . 

The Fourier transform on Li(R d ) and that in L 2 {W l ) are defined as in the 1- 
dimcnsional case, and we denote it by g, if g £ ii(R d ), and by T w [g] , if .9 £ L 2 (R d ). 
The Paley- Wiener functions on R d are given by 

PWl d) := {g £ L 2 (R d ) : F {d) [g] = almost everywhere outside [-7r,7r] rf }. 

The following multidimensional result can be shown in the same way as the 
corresponding result in the 1-dimensional case (see Theorem I2.9|) . Its proof is 
omitted. 
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Theorem 5.1. Let d G N, let X be a fixed postive number, and let [xj : j € N) 

be a Riesz-basis sequence in W 1 . For any f £ PW^ there exists a unique square- 
summable sequence (a(j, A) : j S N) such that 

(57) &(j , l i)jf(i i -ij) = /W, fceN. 

77ie Gaussian Interpolation Operator Jy^ : PWir — > L2[— 7T, 7r] d , defined by 

4 d) (/)(-) = E a ^ A M rf) (-^), 

where (a(j, A) : j G N) satisfies (|57[) . is a well-defined, bounded linear operator from 
PWi d) to L 2 [-ir,ir] d . Moreover, I™ (/) G C (M d ). 

We now generalize Theorems 14.31 and 14.41 to the multidimensional case assuming 
that the underlying Riesz-basis sequence is a grid. For simplicity we restrict our- 
selves to the case when d = 2, but note that our arguments carry over readily to all 
other values of d. It ought to be noted, however, that even this simplistic situation, 
namely that our data sites form a grid, cannot be handled via a straightforward 
multivariate extension of the crucial ingredients from Section 3. In particular, the 
proofs of Corollary 3.7 and Proposition 3.8 do not extend to higher dimensions. So 
we pursue a different tack below. 

We recall some basic tools from tensor products. Let X and Y be two Banach 
spaces and X* and Y* their dual spaces. The algebraic tensor product of X and 
Y is denoted by X®Y, and consists of the vector space of all linear combinations 
of elementary products x®y with x G X and y G Y . All our Banach spaces are 
spaces of functions, and we can therefore identify elements of X®Y with functions 
on a product of sets. 

A norm a(-) on X®Y is called a reasonable cross norm if a(x2y) < ||a;|| • \\y\\, and 
if for <p G X* and V G Y* , the map (f)®ip : X®Y -> C, Y* x i®Vi ^ E <j>(xi)ip(yi), 
is bounded on (X®Y,a), and a{4>®^) = sup u£X0Y . a (u)<i \4>®i>(u)\ < \\<j>\\ ■ If 
a(-) is a norm on X®Y, we denote the completion of X®Y with respect to a(-) 
by X® a Y. 

Proposition 5.2. cf. |Ry| Section and 3.1, Proposition 6.1] 
Let X and Y be Banach spaces. For u G X®Y define 

e{u) := sup{0®7/>(u) : cj> G X*, ip G Y* , \\<f>\\,\\1>\\ < 1} 

(i) e(-) is a reasonable cross norm, and we call e the injective tensor norm on X®Y . 

(ii) If ct(-) is any reasonable cross norm on X x Y , then 

e(u) < a(u), u G X®Y. 

(Hi) For Banach spaces V and W , and bounded operators S : X — > U and T : Y — > 
W, the map S^T : X®Y -> V®W , defined by 9gT(J% =1 x^n) = Y% =1 S(x i yg)T(y i ), 
extends (uniquely, by density) to a bounded operator S(S E T : X® e Y — > V® e W , and 
||S® £ T|| = ||S|H|T||. 

On tensor products of Hilbert spaces we can define a unique tensor norm for 
which the completions are again Hilbert spaces. 
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Proposition 5.3. cf. Kli p. 125 ff] 

Let H and K be Hilbert spaces. 

(i) There is a unique inner product (•, -)h®k on H®K for which 

{hi®ki,h 2 ®k 2 ) H&i = {h\,h 2 ) H {ki,k 2 ) K , h\,h 2 EH and k x ,k 2 EK, 

where (•, -)h and (•, -)k denote the inner products in H and K , respectively. We 
denote the completion of H®K with respect to the corresponding Hilbert norm \\-\\h 

1/2 

(i.e., \\u\\ h = (u,u)^ m , for u E H®K) by H ® h K. 

(ii) || • II/, is a reasonable cross norm on H®K. 

(in) If V and W are two Hilbert spaces, and S : H —> V and T : K : W are two 
bounded operators then the map S(&T : H®K — > V®W , defined via S(®T(Y^ =1 Xi® 
Ui) = X)™=i S( x i)®T(i)i), extends to a bounded operator S®hT : X®hY — > V<S>hW, 
ond||S<8hT|| = ||S||-||T||. 

Let (a : i E N) C H and (/< : i E N) C K. Then (e* : i E N) and (/, : i £ N) are 
Riesz bases, or orthogonal bases for H and K , respectively, if and only if (e,® fj : 
i,j E N) is a Riesz basis, or an orthogonal basis, of H®hK . 

If (CI, S, n) and (f2' , £' , //) are two measure spaces, and H = L 2 (p), and K = 
L 2 (n'), then H®hK is, via the identifcation of f ® g with the function Q x f2' 9 
(u>,u>') i — ► g(u))f(uj l ), isometrically isomorphic to L 2 (/i x //). 

Using Propositon 15.31 and the identification of L 2 (M.) (g> L 2 (M.) with L 2 (M. 2 ), we 
deduce that ^"^'[-J = .?"[•] <8>h.F[-]. Consequently, 

PW^ = (^)- 1 ( J L 2 [- 7 r,^] 2 ) = (^Sfc^-^LahTT.Trl^hLal-Tr.Tr]) 

= PW n ® h PW„. 

For the remainder of this section we fix two Riesz-basis sequences (xj : j E Z) 
and (yj : j E Z), and we put zn t k) = {xi,y k ) for I, k E Z. By Proposition 15. 3[ 
(z(7 : I, k E Z) is a Riesz-basis sequence on M. 2 , and we denote the 1-dimensional 
Gaussian interpolation operator corresponding to A > and (xj : j E Z) by 7 A , and 
the one associated to A > and (yj : j E Z) by /'. The 2-dimensional Gaussian 

(2) 

interpolation operator corresponding to (zn^) '■ l> k E Z) is denoted by 7^ . 
Proposition 5.4. For any positive number X, 

Ix ] =h®hl' x . 

Proof. From Theorem 12.91 Theorem 15. R and Proposition 15. 3[ we find that the 
operators I x ®hl' x and 1^ are bounded linear operators from PW^ into L 2 [— n, ir] 2 . 
So we only need to verify that, 

i ( x\f®h) = i x {f)® h r x {h), f,hEPw^. 

Recall from Theorem 12.91 that there exists a unique pair of square-summable se- 
quences (a(j, A) : j E Z) and (b(j, A) : j E Z) such that 



^2a(j,X)gx{x k - Xj) = f(x k ) and ^b(j,\)g x {y k - Vj) = h(y k ), k E Z. 
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This yields for k, I G Z 

f(x k )h(yi) = Y a(j, \)g\{x k - x 3 ) ^ b(j, \)g x (yi - y s ) 

jez jez 

= X] a (j^) b ( m >^)9\( x k ~ Xj)g\(yi - y m ) 

Therefore the uniqueness in Theorem 15.11 implies that, for every (x,y) €E M 2 , 
I^\f®h){x,y)= a(j,X)b(m,X)g x 2 \(x,y)-(x j ,y m )) 

= a(j,\)b(m,X)gx(x-x j )gx(y-y m ) = I x (f)(x)I' x (h)(y), 

and this finishes the proof. □ 

Proposition 15.41 allows us to transfer Theorems 14.31 and 14.41 to the multidimen- 
sional case. 

Theorem 5.5. Suppose that F e PW^ . Then lim \\l{ 2) (F ) - F\\ Lo(m = 0, and 
lim I^\f)(z) = F(z) uniformly for z £ R 2 . 

(2) 

Proof. Let F G PW^ . The first assertion follows from the aforementioned fact 
that PWir^PWi; is dense in P\v¥\ Proposition [OJ and Theorem fQl 

In order to show the second statement we first note (cf. |Ry| page 50]) that 
the injective tensor product Co(R)<S> £ Co(K) is, via the natural map, isometrically 
isomorphic to Co(K 2 ). By Theorem 14.41 and Proposition [521 (hi) the operators 

I x ®el' x ■ PW«® £ PWv -> C (M 2 ) 

are uniformly bounded. So Propositions l5.2f ii) and l5.3f ii) imply that I x ®h,I[ = I\ 
are also uniformly bounded operators from PW^ to Co(K 2 ). As PW^f^PW^ is 
dense in Pwi d \ we argue as in the proof of the first statement, and conclude from 

(2) (2) 

Theorem SU that P x '(F) converges uniformly to F, for F G PWi > . □ 

6. Further Results on univariate Gaussian interpolation 

In this final section we return to univariate interpolation, in order to discuss 
extensions of some results obtained in Section 2. We begin with a general result 
concerning bi-infinite matrices which appears to be folkloric. We have seen two 
articles which cite a well-known treatise for it, but our search of the latter came 
up emptyhanded. A proof of the said result is indicated in [Ja], but for sake of 
completeness and record, we include a fairly self-contained and expanded rendition 
of this argument here. 

Theorem 6.1. Suppose that (A(j,k))j t keZ is a bi-infinite matrix which, as an 
operator on f 2 (Z), is self adjoint, positive and invertible. 

Assume further that there exist positive constants k and 7 such that \A(j. k)\ < 
Ke — k\ j or ever y p a i r j integers j and k. Then there exist constants k and 7 
such that |^4 _1 (s,i)| < Ke^' s ~*' for every s, t G Z. 
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For the proof of Theorem 16.11 we shall need the following pair of lemmata. 

Lemma 6.2. Let (H, (•, •)) be a Hilbert space, and let A : H — > H be a bounded 
linear, self-adjoint, positive and invertible operator. Let R .= I — Tpnr> where L 
denotes the identity. Then R — R* , (x, Rx) > for every x € H , and \\R\\ < 1. 

Remark 6.3. It is a deep fact (cf. |KR| Proposition 2.4.6]) that if A : H — ► H is a 
bounded linear operator such that (x, Ax) G M for all x E H , then A is self-adjoint. 
Thus a bounded linear operator A : H — > H is self adjoint, positive, and invertible 
if and only if 

inf (x, Ax) > 0. 

x£H,\\x\\ = l 

Proof of Lemma \6. 6 A The symmetry of R is evident. If ||x|| = 1, then 

/ A 

(x,Rx) = \\x\\ 2 - (x, p|-x 

By the assumptions on A and the BCS inequality we see that the term on the right 
of the preceding equation is between and 1. Therefore ||i?|j = sup{(x, Rx) : ||x|| = 
1} < 1. If ||i?|| = 1, then there is a sequence (x n ■ n € N) such that ||a; n || = 1 for 
every n, and 

1 = lira (x n ,Rx n ) = Urn ( 1 - (x n , -rr^Xr 
which contradicts the invertibility of A. □ 

Lemma 6.4. Suppose that (i?(s, £)) s> tez is a bi-infinite matrix satisfying the follow- 
ing condition: there exist positive constants C and 7 such that \R(s,t)\ < Ce -7 ' 15- *' 
for every pair of integers s and t. Given < 7' < 7, there is a constant (7(7,7'), 
depending on 7 and 7', such that \R n (s,t)\ < C n C(-/, 7 / )"~ 1 e~ 7 s ~* for every 

s,te z. 

Proof. Suppose firstly that s ^ t 6 Z, and assume without loss that s < t. Note 
that 

00 t 
(5g) e -7l*-«l e -7'l*-«l _ e -7(«-s) e -T'(t-«) 

s — 1 oc 



+ e-^ s -^e-^' (*"«)+ X! 

U=— OO M=t + 1 

:Si + E 2 + S 3 . 



g-7(u-s)g-7'(«-t) 



Now 



(59) Ei = e -7'(*-«) e (7-7')s^ 



e — u(7-7 7 ) 



g-7'(*-s)g(7-7') s p -s{j-l')-' «(7-7') 
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Moreover, 

e -7'(*-s) 



(60) S 2 =^ e -7" e -7'(t-^) =e -7'(i-s)^ e ^(7+7') < _ , . 

V— 1 V— 1 

whereas 

oo oo c'(t — s) 

(61) E 3 = ^ e -V* e -7(f+t- a ) = e -7(*-s) ^ e -(7+7')« < & _ (7+7<) . 

u=l i; — 1 

If s = t, then 

oo oo ~ 

(62) V e -7k-«l e -7'l*-«l = V" e -(7+7')l«-«l < 

V ' Z-^ ^ - 1 _ e -(7+7') 

li — — oo u— — oo 

From (j58])-(l62l) we conclude that 

l* 2 M)l < c \_ e \^ + t,^) ] = : c9 ^-^- 

The general result follows from this via induction. □ 

Proof of Theorem 15.51 Let R = I — ply be the matrix given in that lemma. As 

r>( ■ A(j,k) \\A\\-A(k,k) 
R{3,k)= lij^k, and R(k,k) = ^ , 

there is some constant C such that \R(s, t)\ < Ce^ 7 ' 5- *' for every pair of integers s 
and t. As A = \\A\\ (I — R), and r := ||i?|| < 1 fLemma l6.2[) . the standard Neumann 
series expansion yields the relations 



(63) A' 1 = \\A\\- 1 ^rR r ' 

n=Q 
N-l 

= \\Al\- 1 Rn + \\A\\- 1 R N Y, Rn = Mir 1 H R" +R N A-\ 



n=0 

N-l oo N-l 



for any positive integer N. As i?°(s, i) = I(s, t) = if s ^ t, s, t G Z, we see from 
(1631) that 



JV-l 



(64) A _1 (M) = || AH" 1 J] iT(s,i) + [R N A' 1 ] (s,t), s^t. 

n=l 

Choose and fix a positive number 7' < 7, and recall from Lemma 16.41 that there is 
a constant (7(7,7') such that \R n (s,t)\ < C"C(7, 7 / )"~ 1 e~ 7 I s- *' for every positive 
integer n, and every pair of integers s and t. So we may assume that there is some 
constant D := 15(7,7') > 1 such that \R n (s,t)\ < D n e~ J for every positive 
integer n, and every pair of integers s and t. Using this bound in (|64p provides the 
following estimate for every s 7^ t: 

N-l 

(65) \A-\s, t)\ < \\A\\- 1 e-' / '^ ^ D n + \\ A' 1 ]]^ 

n=l 

n N 

< \\A\\- 1 e-'' \ s - t \^-j + \\A- 1 \\r N . 
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Let m be a positive integer such that e 7 D 1 /™ < 1, and let s, t £ Z with \s— t\ > 
Writing \s — t\ = Nm + k, < k < m — 1, we find that 



(66) e -i'\»-t\D N = e-t'D^+^m < [ e -f'D 1/m ]^. 

Further, 



s-t 



(67) r 



JV 



| s -t| 



< r 



l/2ml|s-t| 



and combining (|66p and (|67j) with (|65p leads to the following bounds for every 
\s — t\ > m and an appropriately chosen 7 > 0: 

(68) |A- 1 (s,t)| < ML^.[ e -i' D ^/^-t\ + HA- 1 !!^ 1 / 2 " 1 ]! 5 -*! = 0(e-^ s -^). 

On the other hand, if |s — i| < to, we obtain 

(69) < < (||A- 1 ||e m7 )e- 7|s -* 1 , 

and combining ([55)1 with finishes the proof. □ 

A direct consequence of the previous theorem is the following: 

Corollary 6.5. Let A > 0, and let (xj : j 6 Z) be a sequence of real numbers 
satisfying the following condition: there exists a positive number q such that Xj+i — 
x j ^ 9 for every j £ Z. Let A := A\ be a bi-infinite matrix whose entries are given 
by A(j, k) := e~^ Xj ~ Xk ^ , j,k S Z. Then there exist positive constants f3\ and 71, 
depending on A and q, such that |A -1 (s,t)| < Pie"' 1 ^^ , s,t £ Z. 

Proof. The hypothesis on the Xj's implies that \xj —Xk\ > |i — for j, fc G Z. □ 

Remark 6.6. The foregoing result implies, in particular, that A -1 is a bounded 
operator on every £ P (Z), 1 < p < 00. 

We turn now to the second half of this section, in which we introduce the funda- 
mental functions for Gaussian interpolation (at scattered data sites), and set out 
some of their basic properties. 

Theorem 6.7. Let A > be fixed, and let (xj : j £ Z) be a sequence of real numbers 
satisfying the following condition: there exists q > such that Xj+i — Xj > q for 
every integer j. Let A = A\ be the bi-infinite matrix whose entries are given by 
k) = e~ x ( X: >~ Xh ' , j,k £ Z. Given I £ Z, let the l-th. fundamental function be 
defined as follows: 

L hX (x) := L t (x) := ^ A- 1 (fc,/)e- A( ^ a;fc)2 , x £R. 
The following hold: 

(i) The function L\ is continuous throughout M. 

(ii) Each L\ obeys the fundamental interpolatory conditions Li{x m ) = 6 lm , m £ Z. 
(Hi) If in addition, there exists a positive number Q such that Xj+i — Xj < Q for 
every Q, then there exist positive constants P2 and p, depending on X, q, and Q 
such that \Li(x)\ < (32e~ p \ x ~ Xl \ for every x £ R and every I £ Z. 

(iv) Assume that (xj : j £ Z) satisfies the condition stipulated in (Hi). Let (fit : I £ 
Z) be a sequence satisfying the following condition: there exists a positive number 
K and a positive integer P such that \bi\ < for every integer I. Then the 

function K 9 x 1— * X^ez °iLi(x) is continuous on K. 
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Proof, (i) As A^ 1 is a bounded operator on £oo(Z), the sequence (A~ 1 (k, I) : k 6 Z) 
is bounded. Hence the continuity of Li follows from Proposition 12.41 

(ii) Given m € Z, we have 

Li{x m ) = ^2A- 1 {k,l)e- x{Xm ~ Xk)2 = ^2A- l {k,l)A{m 1 k) = I{m 1 l) = 8 lm . 

(iii) The assumption — Xj < Q for every integer j implies that \xk — xi\ < 
\k — l\Q for every pair of integers k and Therefore Corollary 16.51 leads to the 
bound \A~ l {k,l)\ < (3 ie -^ Xk - x '\ , k,leZ, where 72 := ji/Q. Consequently, 

(70) \Li(x)\ <Pi^2 



e -~f2\x k -Xi\ e -\(x-X k ) 



where p := min{A,72}. Therefore 

(71) e p \ x - Xl \\Li{x)\ < /3 1 ^2e pllx ~ x '^ Xk - Xin e- p(x - Xk)2 

fcez 



< f3ij2 ePlx ~ Xkle ~ p(x ~ Xk)2 > xeR - 

fcez 

Fix x € R, and let s be the integer such that x s < x < x s +i. From (j7Tj) we obtain 

s+l 

(72) eH*-*il |£,(a.)| < ft ^ e P\x-x k \ e - P (x-x k ) 2 

k—s 

+ A E 

fcGZ\{s,s+l} 

< 2ft e" Q + ft X! 

tez\{s,s+i} 



e p|:r-:i; fc | e -p(:i;-:t; fc ) 



g p|z-X),l e -p(a:-x*) 



As mg < \x — x s - m \ < (771 + 1)Q for every positive integer m, and (m — l)q < 
\x — x s+m \ < mQ for every positive integer m > 2, the final sum on the right of 
(|72|) is no larger than 

OO OO 

ft ^ eP^+ikg-^V + ft ^ e pm,3 e" p(m ~ 1)2<i2 =: ft. 

m— 1 m— 2 

Combining this with (j?2"|) and (|71 [) finishes the proof. 

(iv) Each summand is continuous by assertion (i). Let and let x s < x < 

x s+ i. Assertion (iii) implies that 



s+l 



Y / b l L l (x)\<Kf3 2 [j2\l\ P + E 

zez i=s iez\{s,s+i} 



\ p e -p\x-x t \ 



s+l 



-P g -P(™-1)9 



s m— 1 



It follows that the series Yliei.biLi(x) is locally uniformly convergent, whence the 
stated result follows. □ 
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The counterpart of Part (iii) of the foregoing theorem for Gaussian cardinal 
interpolation was obtained in |RSlj . and its analogue for spline interpolation was 
proved in lrlalRI 



Earlier in this paper we have discussed Gaussian interpolation operators asso- 
ciated to functions, specifically to bandlimited functions. Here our perspective 
changes slightly, as we begin to think of these interpolation operators acting on 
sequence spaces. 

Theorem 6.8. Let X > be fixed. Suppose that (xj : j G Z) is a real sequence 
satisfying the following condition: there exist positive numbers q and Q such that 
q < Xj+i — xj < Q for every integer j. Let A = A\ be the bi-infinite matrix 
whose entries are given by A(j,k) = e - A (^-z fc ) ^ - k e z ^ Given 

p G [1, oo], and 

y := (yi : I € Z) G £ p (%), define 

I x {y,x) := J2( A ^yh^ Mx ' Xk) \ x e M, 
fcez 

where (A~ l7 y)k denotes the k-th component of the sequence A^ 1 ]]. The following 
hold: 

(i) The function l9a;n L x (y,x) is continuous on R. 

(ii) Lf x is any real number, then 



zez 

where (Li : I G Z) is the sequence of fundamental functions introduced in the pre- 
ceding theorem. 

(iii) There is a constant fa, depending on X, q, Q, and p, such that 

ll^(y,0IUp(R) ^ A5||y|k(z), 

for every y G £ P (Z). 

Proof, (i) As A^ 1 is a bounded operator on £ p (Z), the sequence A~ x y is bounded. 
Hence the continuity of l x (y, ■ ) follows from Proposition 12.41 
(ii) Let Then 



e -\(x-x k ) 2 



fcez k& 



The required result obtains by interchanging the order of summation, which is 
justified by the following series of estimates, the first of which is consequent upon 
Corollary 16. 51 and the last of which follows from Proposition [ 



^[^M- 1 (fe I 0l]e- A ^ )2 =o(^[^e^l' 

fcez zez fcez zez 

= 0(j2e~ Hx ~ Xk)2 ) =0(1) 



e -\(x-x k y 
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(iii) We begin with p — oo. If x s < x < a; s +i, assertion (ii), the triangle 
inequality, and a now familiar argument lead to the following: 



(73) \h(V^)\ 



E^ Lj ( a 



< (3 2 [\y s \ + \y s+1 \ + £ \y s - m \e- pmq + £ \y s+m \e-^ m -^ 
2fo 



< 



1 - e-Pi 



m—l 



ra=1 



It is immediate that \\I x (y, OIU^CK) < #3 1 1 2/1!^ (z)- 

Suppose now that 1 < p < oo. Let J s := [x s , x s + i), s € Z, and recall from ([73 
that, if x S J s , then 



\i x (y,x)\ < fo\\y s \ + \y s+ i\ + E \ys- m \e- pmq +J2\y 



-p(m— l)q 



m=2 



Therefore, as \x s +i — x s \ < Q, we have 
(74) J \I X (y, x)\p dx < QPl [\y s \ + \y s+1 \ 

OO 

+ E 1^- 1 



-pmq 



El*- 



-p(m-l)g 



Let u = (uk ■ k € Z) and i> = : fc € Z) be a pair of sequences defined by 
w/c = \Vk\, k e Z, and 



I if G {—1,0}; 

„ fc = ^ e -pfe<7 if fc > 1; 
e p(fc+i)« if k < _ 2 



Then (|74|) may be recast as follows: 

\I x (V,x)\*dx<Q% 



m£Z 



l s—m u m 



so that 
(75) 



I^a(2/»-)||l p (i;) 



]T / \h{y,x)\ p dx 



1 1/P 



< 1/p /3 2 



E 



E 



Now the Generalized Minkowski Inequality (cf. (HLP[ p. 123]) implies that 



(76) 



E 

.sGZ 



E 



l s—ra u m 



p-i 1/p 



< 



E 



Ei 



i/p 



< 



x _ e - pq llylU P (z)' 



and a combination of ([76) and ([75]) completes the proof. 



□ 
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We conclude the paper with a few remarks. Suppose that / is a bandlimited 
function, and let dk = f(xk), k € Z. We have seen in Proposition 12.81 that this 
sequence d = (dj~ : k E Z) is square summable. Furthermore, as observed in the 
course of the proof of Theorem I2.9f v). the Gaussian interpolant I x (/)(•) studied 
earlier coincides with the function I x (d, •) introduced in this section. Thus the 
final part of the previous theorem presents a twofold generalization of the estimate 
(f2"5| : to values of p E [1, oo] other than 2, whilst requiring only that the underlying 
set of sampling points satisfies condition (jTUJ) . In particular we do not assume in 
Theorem I6.8f iii) that (xj : j € Z) is a Riesz-basis sequence. However, it is not 
without interest to note that the main convergence theorems obtained in Section 4 
do not hold for data sites which merely satisfy the quasi- uniformity condition (|10|> . 
For example, let X := Z \ {0}, and let / be the bandlimited function given by 
f(x) := sin(7T2:)/(7ra), x E R. As f(xj) — for every Xj € X, I x (f) is identically 
zero, so it is manifest that I x (f) does not converge to /. 

Counterparts of the result obtained in part (iii) of the preceding theorem, for 
the case when xj — j, may be found in RS2J and RSI . However, those estimates 
are much more precise in nature. 

Suppose that (xj : j € Z) is a strictly increasing sequence of real numbers 
satisfying the following condition: 

(77) \ Xj - j\ < D < 1/4, jeZ. 



Then (xj) is a Riesz-basis sequence [Kaj . Let (Li t \ : I £ Z) be the associated 
sequence of fundamental functions defined in Theorem 16. 71 Define 



and let 



G(x) := (x - s ) I] - f: J ^ - ) . 

G ^--= (x-%(x iy X ^ l ^ 
It is shown in ILej each Gi is a bandlimited function satisfying the interpolatory 
conditions Gi(x m ) — S lm , m 6 Z. Thus we find that I x (Gi) = L/.a, and conclude 
from Theorems 14.31 and 14.41 that lim A _ +0 + Li : \ = Gi in Z^QR) and uniformly on R. 
As a particular example, we learn from |Hi| that, if xg =0 and Xj = X-j — j + c 2 /j, 
j > 1, < |c| < 1/2, then (xj) fulfills (1771) . and that the corresponding function G 
is given in closed form: 

G(x) = x[cos(ir(x 2 - 4c 2 ) 1/2 ) - cos(7ra;)]/2 sinh(Trc), x E R. 
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